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Abstract 

The high-pressure structural behaviour of the flouroperovskite KMgF3 is investigated by the- 
ory and experiment. Density functional calculations were performed within the local density ap- 
proximation and the generalized gradient approximation for exchange and correlation effects, as 
implemented within the full-potential linear muffin-tin orbital method. In situ high-pressure pow- 
der x-ray diffraction experiments were performed up to a maximum pressure of 40 GPa using 
synchrotron radiation. We find that the cubic Pm3m crystal symmetry persists throughout the 
pressure range studied. The calculated ground state properties - the equilibrium lattice constant, 
bulk modulus and elastic constants - are in good agreement with experimental results. By ana- 
lyzing the ratio between the bulk and shear modulii, we conclude that KMgF3 is brittle in nature. 
Under ambient conditions, KMgF3 is found to be an indirect gap insulator with the gap increasing 
under pressure. 



I. INTRODUCTION 



KMgF3 is a technologically important flouroperovskite. For example, it is used as a 
vacuum-ultraviolet-transparent (VUV-transparent) material for lenses in optical lithography 
steppers 1 and in electro-optical applications^. When doped with lanthanide ions, it is 
a very promising material for scintillators^ and radiation dosimeters^. In addition the 
physical properties of KMgF 3 may have implications for understanding of the Earth's lower 
mantle^. 

KMgF 3 was first synthesized by van Arkel^ and has a simple cubic perovskite structure at 
room temperature^. KMgF 3 demonstrates great stability under high compression and has 
not been found to undergo any phase transition at any temperature or pressure, suggesting 
it may be used as an internal X-ray calibrant 11 . 

Several experimental studies of the ground state properties of KMgF 3 have been per- 
formed. The elastic constants at ambient pressure have been measured by Rosenberg and 
Wigmore 12 and by Reshchikova,— while Jones investigated their pressure and temperature 
dependence^. 

From the theoretical side, electronic structure calculations for KMgF 3 have been carried 
out by means of linear combination of atomic orbitals 15 , including the effects of doping with 
transition metal impurities in KMgF 3 16 . The electronic structures of divalent 3d transition 
metal impurities doped in KMgF 3 have been investigated by the pseudopotential method^, 
and the properties of vacancies were studied by Hartree-Fock cluster calculations 18 . The 
structural, electronic and optical properties of KMgF 3 were recently investigated by the 
full-potential linear augmented plane wave (FP-LAPW) method 1 -. 

The present work is a combined theoretical and experimental study of the ground state 
and high-pressure properties of KMgF 3 . We present the equation of state resulting from 
high-pressure diamond-anvil cell experiments on KMgF 3 up to 40 GPa. We also present 
the equation of state, the elastic constants and the electronic structure from theoretical 
calculations using two different approximations for the exchange-correlation functional. 

The remainder of the paper is organized as follows. Details of the computational method 
as well as details of the experimental setup are outlined in section 2. The measured and 
calculated equations of state are presented in section 3 together with calculated ground state 
properties and elastic properties. The electronic structure and the pressure variation of the 
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band gap are discussed in section 4. Finally, conclusions are given in section 5. 



II. COMPUTATIONAL AND EXPERIMENTAL DETAILS 



A. The electronic structure method 

The all-electron full-potential linear muffin tin orbital (FP-LMTO) method 2 ^ is used to 
calculate the total energies and basic ground state properties of KMgF 3 presented here. 
In this method, the crystal volume is split into two regions: non-overlapping muffin-tin 
spheres surrounding each atom and the interstitial region between the spheres. We used a 
double k spdf LMTO basis (each radial function within the spheres is matched to a Hankel 
function in the interstitial region) to describe the valence bands. In the calculations we 
included the 3s, 3p, 4s, 4p, and 3d bases for potassium, the 3s, 2p, 3p, and 3d bases for 
magnesium, and the 2s and 2p bases for fluorine. The exchange correlation potential was 
calculated within the local density approximation (LDA)^ as well as the generalized gradient 
approximation (GGA) scheme 2 ^. The charge density and potential inside the muffin-tin 
spheres were expanded in terms of spherical harmonics up to / maa; = 6, while in the interstitial 
region, they were expanded in plane waves, with 14146 waves (energy up to 156.30 Ry) 
included in the calculation. Total energies were calculated as a function of volume for a (16 
16 16) k-mesh containing 165 k-points in the irreducible wedge of the Brillouin zone and 
were fitted to the Birch equation of state 23 to obtain the ground state properties. 



The elastic constants were obtained 
conserving strains, as outlined in Refs. 



rom the variation of the total energy under volume- 
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B. Experimental details 

The high-pressure x-ray diffraction measurements used a sample of polycrystalline KMgF 3 
prepared by the solid state reaction method from high purity constituent materials as de- 
scribed elsewhere in several reports 2 ^ 2 ^ 1 ^. Diffraction patterns collected at ambient temper- 
ature and pressure showed a cubic (Pm3m) symmetry with a cell parameter a = 4.0060(2) 
A for KMgF 3 which closely agrees with earlier reports 3 ^*. High pressures were generated by 
a Merrill-Bassett type diamond-anvil cell (DAC). An 185-/xm sample chamber was formed 
in a rhenium metal gasket with a pre-indention of 60-/mi thickness. The powder sample 
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FIG. 1: Powder x-ray diffraction patterns recorded at various pressures up to 32 GPa. The 
indexing in terms of the simple cubic structure is given. 



was loaded in the gasket with a few ruby grains and silicone fluid as pressure transmitting 
medium 3 -^. Diffraction experiments were performed at the 16ID-B undulator beam line of 
the High Pressure Collaborative Team (HPCAT) of the Advanced Photon Source (APS). A 
monochromatic x-ray beam with a wavelength of 0.4218 A was focused down to a size of 
30 x 30 /im. Diffraction images were collected with an image plate detector for an exposure 
time of 10 sec. The distance between the sample and the detector and the inclination angle 
of the image plate were calibrated using a Ce02 standard. 

The two dimensional images were subsequently integrated to one dimensional diffraction 
patterns using the Fit2D software^ 2 -. The cell parameters were obtained by analyzing the 
diffraction patterns with the JADE software package and the P-V data obtained was fitted 
with a second-order Birch-Murnaghan 23 equation of state. The standard ruby fluorescence 
technique and the newly proposed ruby pressure scale of Holzapfel 33 were used to obtain the 
pressures in the sample chamber. 
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III. GROUND STATE AND ELASTIC PROPERTIES 



Powder x-ray diffraction patterns collected at several pressures are shown in Figure 1. On 
compression, the diffraction patterns remain unchanged up to 40 GPa, except for the shifts 
of diffraction lines caused by the decreasing lattice constant. This implies that no structural 
transformations occur up to 40 GPa in KMgF3. Figure 2 shows the measured equation of 
state of KMgF 3 and compares it with theoretical curves calculated within the LDA and GGA. 
The better theoretical description is obtained with the LDA which is somewhat surprising, 
since usually the GGA provides an improvement over LDA. At low pressures the LDA volume 
is slightly smaller than the experimental one, while the situation reverses at high pressures, i. 
e., altogether the LDA predicts KMgF 3 to be stiffer than experimental observations. GGA 
on the other hand greatly overestimates the equilibrium volume at ambient pressure, which 
is the main reason for the poor agreement with experiment. If the GGA curve is scaled 
througout the pressure range with the error in equilibrium volume at P = 0, nearly perfect 
agreement is found with experiment (not shown). 

The lattice constant and bulk modulus measured in the present work as well as values 
calculated within the LDA and GGA approximations are given in Table I. Results from 
earlier experimental and theoretical works are quoted for comparison. The bulk modulus 
obtained in our experiments B = 71.2(2) GPa with B' = 4.7(3) compares well with other 
experimental results listed in Table 1 and also with NaMgF 3 reported recently by Liu et 
al^I (_B = 76.0(1.1) GPa). The lattice constant obtained within the LDA is 1.1 % lower 
than the experimental value, while the corresponding bulk modulus is 22% higher than the 
experimental value, which is the usual kind of accuracy of LDA. However the calculated 
LDA lattice constant from the present work agrees quite well with the experimental work 
when compared to the earlier FP-LAPW(LDA) calculations in which the reported lattice 
constant is 2.4% lower than the experimental valued. The LDA bulk modulus obtained 
from the present calculation agrees quite well the published FP-LAPW(LDA) results. When 
comparing the results obtained within GGA, the lattice constant is 1.9 % higher than the 
experimental value, whereas our results for the bulk modulus is within the spread of the 
experimental data. This truly excellent agreement regarding the bulk modulus is, however, 
a bit fortuitous. Since the calculated equilibrium volume is overestimated with GGA (and 
underestimated with LDA), an error - solely depending on the error in volume - is introduced 



in the calculated bulk modulus. Therefore, we recalculated the bulk modulus also at the 
experimental volume in a manner similar to our earlier work^ (see Table I). We find that this 
diminishes the discrepancies between the LDA and GGA results, as expected. In addition, 
the LDA bulk modulus now becomes smaller than the GGA one for KMgF 3 , and both 
functional are seen to actually overestimate the bulk modulus, by approximately 14 % 
(LDA) and 34 % (GGA). 

The present experiments on KMgF3 relate to recent experiments performed for NaMgF3 
and alloys of NaMgF 3 and KMgF 3 . The crystal chemistry of Nai_, r K r MgF 3 and NaMgF 3 
was studied in detail at ambient and at high pressures by Zhao et al.— >2*. NaMgF 3 under- 
goes a reversible phase transition from orthorhombic (Pbnm) to tetragonal (P4/mbm) and 
then to cubic structure (Pm3m) upon compression. These phase transitions require either 
compositional changes, by increasing the K concentration to 40%, or changing temperature 
or pressure. The structural changes in these perovskites are due to octahedral tiltings and 
shortening of Mg-F bonds compared to the cubic phase. A direct transformation from or- 
thorhombic to cubic structure in NaMgFs, however, requires a very high temperature (1038 
K) . Moreover, the transition temperature is reported to increase with pressure. The temper- 
ature dependence of the crystal structure of KMgF 3 was recently investigated by neutron 
powder diffraction by Wood et al.— from 4.2 K to 1223 K, and the cubic symmetry was 
found to be stable throughout this temperature range. The thermal expansion as well as 
the atomic displacement parameters obtained in their experiments show that the F ions 
behave less anistropically than in NaMgF 3 at such high temperatures. On comparing these 
results with the present high-pressure diffraction experiments on KMgF 3 , one may speculate 
that application of pressure alone would not be sufficient to induce structural changes in 
KMgF 3 as the cubic phase is very stable. Such a phase transformation if any, would require 
either application of very high temperature or a composition change in the system to achieve 
changes in the order parameters. Asbrink et al.— have studied single crystals of the transi- 
tion metal bearing perovskite KMnF 3 , which is isostrucural to KMgF 3 , under high pressure 
and observed a cubic-to-tetragonal phase transition at a critical pressure of 3.1 GPa. On 
combining these results, phase transitions from the cubic symmetry may be expected with 
a combination of composition change, temperature and pressure in KMgF 3 . A systematic 
study on the octahedral tilting and order parameters with other dopant compositions and the 
effect of external thermodynamical variables are further required to understand the phase 
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FIG. 2: (Color Online) Equation of state of KMgF3 in the pressure range from 0-40 GPa. The ex- 
perimental datapoints are marked by circles, while theoretical results obtained with LDA and GGA 
are shown as full (blue) and dashed (red) curves. Volumes are given relative to the experimental 
equilibrium volume Vb = 64.288 A. 

stability of KMgF 3 . 

The elastic constants of KMgF3 calculated within LDA and GGA are listed in Table 2 
where they are also compared to experimental results as well as earlier calculations. The 
LDA overestimates all of the Cn, Cyi and C44 elastic constants by between 10% and 22% 
compared to experiment.— The elastic constants obtained within GGA are much closer 
to the experimental values than are the LDA results. For instance, both C\\ and C12 are 
within the experimental spread. Of course, the elastic constants also depend sensitively on 
the volume, and, therefore, the same argument as for the bulk modulus can be applied here. 
We have, however, refrained from recalculating all the elastic constants with the volume 
correction, but wish to mention that the excellent agreement between experiment and the 
GGA elastic constants should be interpreted with care. Another point of caution is the fact 
that the calculated values pertain to kelvin, while experiments are performed at room 
temperature. Finite temperature generally tends to reduce the elastic constants because 
of thermal expansion. Using the calculated elastic constants we calculated the anisotropy 
factor A = 2C M /(C n - C 12 ). We find an A = 0.91 for LDA and A = 1.12 for GGA. 
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The experimental result is 1.05, measured at room temperature^ 1 ^ 1 ^, which is closer to 
but lower than the GGA value. However, the anisotropy factor is found to decrease as the 
temperature is lowered^. 

A simple relationship, which empirically links the plastic properties of materials with their 
elastic moduli was proposed by Pugh^. The shear modulus G represents the resistance to 
plastic deformation, while the bulk modulus B represents the resistance to fracture. A high 
B/G ratio is associated with ductility whereas a low value corresponds to brittle nature. The 
critical value which separates ductile and brittle materials is around 1.75, i.e. if B/G > 1.75 
the material behaves in a ductile manner, otherwise the material behaves in a brittle manner. 
Frantsevich 42 , in a similar fashion has suggested B/G ~ 2.67 as the critical value separating 
brittle and ductile behavior. In the case of KMgF 3 the calculated value of B/G is 1.5 within 
LDA and 1.4 within GGA, hence classifying this material as brittle. 

Pettifor— suggested that the angular character of atomic bonding in metals and com- 
pounds, which also relates to the ductility, could be described by the Cauchy pressure 
C\2 — C44. For metallic bonding the Cauchy pressure is typically positive. On the other 
hand, for directional bonding with angular character, the Cauchy pressure is negative, with 
larger negative pressure representing a more directional character. These correlations have 
been verified for ductile materials such as Ni and Al that have typical metallic bonding, as 
well as for brittle semiconductors such as Si with directional bonding^. In the ionic com- 
pound KMgF3, the calculated Cauchy pressure is -10 GPa within LDA and -15 GPa within 
GGA, in good agreement with the nonmetallic characteristics of KMgF3. 

Table 3 presents sound velocities as derived from the calculated elastic constants^. The 
calculated sound velocities agree quite well with the experiments, in particular for the GGA 
values, which is a consequence of the somewhat fortuitous good agreement between the 
measured and GGA calculated elastic constants. 

IV. ELECTRONIC STRUCTURE 

The calculated electron band structure of KMgF 3 is shown in Figure 3 with the ensuing 
density of states in Figure 4. The valence bands consist of the F p bands with a gap of 
7.24 eV to the conduction band, which is dominated by K states. The LDA bands are 
almost identical, however with a gap of only 6.95 eV. The gap increases almost linearly with 
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FIG. 3: Band structure of KMgF3 (using GGA, at the experimental lattice constant). The zero of 
energy is set at the position of the conduction band minimum. 



compression, at the rate 



v dE, 
dV 



-7.1 eV. 



The conduction band minimum occurs at the T point, while the valence band maximum 
occurs at the R-point (1/2, 1/2, 1/2)—. The largest occupied energy level at the M-point 
(1/2,1/2,0)^ is marginally lower (by ~ 0.02 eV) than the valence band maximum at R, 
and it remains lower throughout the pressure range studied here. 

V. CONCLUSIONS 

In conclusion, we find that pure cubic KMgF3 is very stable under high compression. 
From our analysis we also find that it is a brittle system and an indirect gap insulator whose 
gap increases with pressure 
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TABLE I: Calculated lattice constants (in A), Bulk modulus Bq (in GPa) and its pressure derivative 
B' , of KMgF3 at the theoretical equilibrium volume compared with the experiment and other 
theoretical calculations. The bulk moduli have been calculated both at the experimental and 
theoretical volume (Bq(Vq XV> ) and Bo(VQ h ), respectively) 

Lattice constant B (Vj h ) B (Vq XP ) B' 



GGA a 4.0809 72.01 97.85 4.65 

LDA a 3.9630 91.47 83.23 4.79 

LDA, LAPW-? 3.91 90.97 - 4.64 

Expt. 4.0060(2) a , 3.973 6 , 3.978±0.05 c , 71.2(2) a ,70.4^, 4.7(3) a 

3.9897 d , 3.993 e , 3.9839 / 75.1 / \75.6 i 
a Present work, 6 Ref. 10 , r R<4. :il . d ReL m , 'RclV'"'. ^Ref, 29 , ^Ref, 12 , ^Ref, 13 , *Ref. 14 iReiM 

TABLE II: Calculated elastic constants, shear modulus (G), and Young's modulus (E) all expressed 
in GPa, and Poisson's ratio v of KMgF3 at the theoretical equilibrium volume 
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TABLE III: Calculated longitudinal, shear, and average wave velocity (vi, v s , and v m , respectively) 
in m/s for KMgF3 at the theoretical equilibrium volume 

Vi v s v m 



Present 



Expt. 



LDA 
GGA. 



7402 
6706 
6470 a ,6540 6 



4396 
4073 
3940 a ,3900 b 



4870 
4507 
4290 c 



4^U 



a : Wave vector along < 100 > direction, Refll2l. : Wave vector along < 110 > direction, Refll2l. 

c : Ref 
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